Method for determination of optimum fields of permanent magnet structures with linear magnetic characteristics

ABSTRACT

The invention is directed to a method for determining the fields of permanent magnet structures with a surface or boundary solution method for the magnetic material with linear characteristics with small susceptibility and large permeabilities of the ferromagnetic materials.

FIELD OF THE INVENTION

This invention relates to an improved method for determining the optimum fields of permanent magnetic structures having linear magnetic characteristics, for enabling the more economical production of magnetic structures.

BACKGROUND OF THE INVENTION

Exact solutions can be achieved in the mathematical analysis of structures of permanent magnets under ideal conditions of linear demagnetization characteristics and for some special geometries and distributions of -magnetization. For instance, an exact mathematical procedure can be followed to design a magnet to generate a uniform field in an arbitrarily assigned polyhedral cavity with perfectly rigid magnetic materials and ideal ferromagnetic materials of infinite permeability.

In general, for arbitrary geometries and real characteristics of magnetic materials, only approximate numerical methods can be used to compute the field generated by a permanent magnet. The capability of handling systems of a large number of equations with modern computers has led to the development of powerful numerical tools such as the finite element methods, in which the domain of integration is divided in a large number of cells. By selecting a sufficiently small cell size, the variation of the field within each cell can be reduced to any desired level. Thus the integration of the Laplace's equation in each cell can be reduced to the dominant terms of a power series expansion and the constants of integration are determined by the boundary conditions at the interfaces between the cells. An iteration procedure is usually followed to solve the system of equations of the boundary conditions and the number of iterations depends on the required numerical precision of the result.

In applications where the field within the region of interest must be determined with extremely high precision, the large number of iterations may become a limiting factor in the use of these numerical methods. It is beyond the scope of this disclosure to provide a detailed explanation of past techniques for this purpose.

A special situation is encountered in magnetic structures that make use of the rare earth permanent magnets that exhibit quasi linear demagnetization characteristics with values of the magnetic susceptibility small compared to unity. A magnetic structure composed of these materials and ferromagnetic media of high magnetic permeability can be analyzed with a mathematical procedure based on a perturbation of the solution obtained in the limit of zero susceptibility and infinite permeability.

Structures composed of ideal materials of linear magnetic characteristics present a special situation where an exact solution is formulated by computing the field generated by volume and surface charges induced by the distribution of magnetization at the boundaries or interfaces between the different materials.

SUMMARY OF THE INVENTION

The determination of the field in this ideal limit can be developed with a boundary solution method which may be formulated in a way that substantially reduces the number of variables as compared to the finite element method. The invention is therefore directed to a method for determining the fields of permanent magnet structures with a surface or boundary solution method for the magnetic material with linear characteristics with small susceptibility and large permeabilities of the ferromagnetic materials.

BRIEF DESCRIPTION OF THE DRAWING

In order that the invention may be more clearly understood, it will now be disclosed in greater detail with reference to the accompanying drawing, wherein:

FIG. 1 illustrates the magnetic conditions at the interfaces of three media;

FIG. 2 defines the most general configuration of the magnetic media;

FIG. 3 illustrates one of the surfaces of FIG. 2;

FIG. 4 illustrates a strip of infinite permeability in a uniform magnetic field;

FIG. 5 is a table showing the distribution of surface charges along the strip for n=20;

FIG. 6 show a plot of equipotential lines generated by the strip;

FIG. 7 shows the equipotential lines when the angle α=O;

FIG. 8 shows the equipotential lines around the strip the angle α=45°;

FIG. 9 illustrates an equilateral hexadecagon at 45° with respect to a uniform field. In this figure the magnetic permeability of the material is infinite;

FIG. 10 illustrates the polyhedron of FIG. 9 assuming μ₀ /μ=0.5;

FIG. 11 illustrates a structure of uniformly magnetized material and zero-thickness plates;

FIG. 12 illustrates the field configuration of the structure of FIG. 11;

FIG. 13 illustrates the field configuration corresponding to the separation of inclined sides;

FIG. 14 illustrates the field configuration within the structure under the condition Φ₃ =Φ₄ =0;

FIG. 15 illustrates the field configuration outside of the structure under the condition Φ₃ =Φ₄ =0; and

FIGS. 16-18 constitute a flow diagram of the method of the invention.

DETAILED DISCLOSURE OF THE INVENTION Field of structure for ideal materials with susceptibility vm=0 and μ=∞.

Consider the structure of FIG. 1 composed of three media: a nonmagnetic medium in region V₁, an ideal magnetic medium of zero magnetic susceptibility (χ_(m) =0) in region V₂, and an ideal ferromagnetic medium of infinite magnetic permeability μ in region V₃. This figure represents the most general interface and defines a basic boundary condition.

Because of the assumption μ=∞, the region V₃ is equipotential and so are the interfaces S₁, S₂ between the region V₃ and the two regions V₁ and V₂. Thus, at each point of interfaces S₁, S₂ the intensities H₁, H₂ of the magnetic field computed in regions V₁ and V₂ are perpendicular to the interfaces, as indicated in FIG. 1.

Assume a unit vector n perpendicular to the boundary surface of region V₃ and oriented outward with respect to V₃. The intensity of the magnetic field induces a surface charge σ on interfaces S₁, S₂ given by

    σ=μ.sub.0 H·n                            (1)

On the interface S₃ between the region V₁ of nonmagnetic material and the region V₂ of magnetic medium, the surface charge density σ₃ is given by

    σ.sub.3 =μ.sub.0 (H.sub.2 -H.sub.1)·n.sub.3(2)

where the unit vector n₃ is perpendicular to S₃ and oriented from region V₁ to region V₂. The magnetic induction B₁ in the region V₁ is

    B.sub.1 =μ.sub.0 H.sub.1                                (3)

and the magnetic induction B₂ in the region V₂ of zero magnetic susceptibility is

    B.sub.2 =J+μ.sub.0 H.sub.2                              (4)

where J is the remanence of region V₂. On interface S₃ vectors B₁, B₂ satisfy the condition

    (B.sub.2 -B.sub.1)·n.sub.3 =0                     (5)

Thus eq. (2) reduces to

    σ.sub.3 =-J·n.sub.3                         (6)

In general, a singularity of the intensity H occurs at the intersection P of the interfaces unless the geometry of the interfaces and the surface charge densities satisfy the condition

    Σσ.sub.h τ.sub.h =0                        (7)

where h are integers and τ_(h) are the unit vectors tangent to the interfaces at point P and oriented in the direction pointing away from the interfaces.

Assume a number N of surfaces S_(h) of μ=∞ media as shown in FIG. 2. This figure illustrates the most general configuration with arbitrary distribution of remanence J. The region is limited by plural regions S enclosing media of given μ. The boundary S₀ limits the region of interest. FIG. 3 illustrates an arbitrary one of the surfaces of FIG. 2, in greater detail. The external region surrounding the N surfaces is a medium of zero magnetic susceptibility with an arbitrary distribution of remanences J, which is equivalent to a volume charge density

    v=-V·J                                            (8)

In the particular case of a uniform magnetization of the external region, the vector J is solenoidal and the distribution of magnetization reduces to surface charges σ_(i) on the interfaces between the regions of remanemces J_(i-1) and J_(i)

    σ.sub.i =(J.sub.i-1 -J.sub.i)·n.sub.i       (9)

where n_(i) is the unit vector perpendicular to the interface and oriented from the region of remanence J_(i-1) to the region of remanence J_(i). Eq. (7) is a particular case of eq. (9).

At each point P of the structure of FIG. 2 the scalar magnetostatic potential is ##EQU1## where V is the volume of the external region, σ_(i) is the surface charge density induced by J at a point of S_(i), σ is the distance of point P from a point of volume V, and σ_(i) is the distance of P from a point of surface S_(i). In the limit μ=∞ the surface charge densities σ_(i) in eq. (10) are determined by the boundary conditions

    Φ(P.sub.h)=Φ.sub.h                                 (11)

where P_(h) is a point of surface S_(h) and Φ_(h) is the potential of surface S_(h). Equation (11) is an identity that must be satisfied at all points of S_(h).

Equations of the type of equations (10) and (11) may be employed in the determination of the magnetic fields of permanent magnetic structures, using a volumetric analysis. This approach, however requires extensive calculations, especially when complex structures are to be analyzed. In accordance with the present invention, as will now be discussed, much simpler and less time consuming calculations may be made employing surface analysis, to thereby reduce the effort required for the production of a magnetic structure having desired characteristics.

By definition, each surface S_(h) immersed in the magnetic field generated by J cannot acquire a non zero magnetic charge. Thus the distribution of surface charges σ on each surface S_(h) must satisfy the condition ##EQU2## Thus, by virtue of eqs. (10) and (11), the unknown quantities σ_(i), Φ_(h) are the solution of the system of equations (12) and the identities ##EQU3## where ρ_(h) is the distance of a point P of surface S_(h) from a point of volume V, and ρ_(h),i is the distance of P from a point of surface S_(i). For i=h , ρ_(h),i is the distance between two points of surface Sh_(h).

In eq. (13) the independent variables Φ_(h) are the potentials of surfaces S_(h) relative to a common arbitrary potential of a surface S₀ that encloses the structure of FIG. 2. In particular S₀ may be located at infinity.

In eq. (13) ρ_(i),h is zero for the element of charge located at the point where the scalar potential is computed. However, as long as σ_(i) is finite, the integral of the left-hand side of eq. (13) does not exhibit a singularity. Consider a circle of small radius r on surface S_(i) with the center at a point P. For r→0, the contribution of the surface charge σ_(i) within the circle of radius r to the potential at P is ##EQU4##

Magnets with Linear Characteristics of Magnetic Media and Arbitrary Ferromagnetic Materials

Eqs. (12) and (13) are based on the assumption of ideal materials characterized by χ_(m) =0 and μ=∞. Assume now that the magnetic material has a linear demagnetization characteristic with a non zero value of the magnetic susceptibility χ_(m)

    χ.sub.m <<1                                            (15)

Assume also a linear characteristic of the ferromagnetic material with a magnetic permeability such that ##EQU5## The magnetic induction in the region of the magnetized material is

    B=J+μ.sub.0 (1+χ.sub.m)H                            (17)

The solution of the field equation within the magnetized material can be written in the form

    B=B.sub.0 +δB, H=H.sub.0 +δH                   (18)

where B₀, H₀ are the magnetic induction and the intensity of the magnetic field in the limit χ_(m) =0. By virtue of eq. (15) one can assume

    |δB|<<|B|, |δH|<<|H.sub.0 |(19)

By neglecting higher order terms, eq. (17) yields

    δB=μ.sub.0 χ.sub.m H.sub.0 +μ.sub.0 δH (20)

i.e., δB and δH are related to each other as if the magnetic material was perfectly transparent (χ_(m) =0) and magnetized with a remanence

    δJ=μ.sub.0 χ.sub.m H.sub.0                    (21)

Thus, the first order perturbation δΦ of the scalar potential is a solution of the equation

    δ.sup.2 (δΦ)=-δ·(χ.sub.m H.sub.0)=-H.sub.0 ·δχ.sub.m -χ.sub.m δ·J                                        (22)

Assume that the magnetic structure is limited by surfaces S_(h) of infinite magnetic permeability materials. By virtue of eqs. (13) and (22), the first order perturbation δΦ and δσ_(h) of the potential and surface charge density on these surfaces are the solution of the identities ##EQU6## and the equations ##EQU7## In the limit 16, the finite magnetic permeability of the ferromagnetic materials inside surfaces S_(h) results in an additional perturbation of the potential in the magnetic structure and in a non zero magnetic field inside surfaces S_(h). At each point of S_(h), the intensities H_(e) and H_(i) of the magnetic field outside and inside the ferromagnetic material satisfy the boundary condition ##EQU8## where n is a unit vector perpendicular to S_(h) and oriented outwards with respect to the ferromagnetic material.

H_(e) and H_(i) are the intensities at two points P_(e), P_(i) at an infinitismal distance from P within the regions outside and inside S_(h) respectively.

The boundary conditions on surface S_(h) will be satisfied by replacing the medium of permeability μ with a surface charge distribution σ on S_(h) and by assuming that:

    μ=μ.sub.0                                            (26)

everywhere. At points P_(e), P_(i) the intensity generated by an element of charge σdσ at P is perpendicular to S_(h) and is given by: ##EQU9## at P_(e) and P_(i) respectively. Thus the normal components of H_(e), H_(i) suffer a discontinuity at P given by: ##EQU10## and because of equation 55 the charge σ(P) satisfies the equation ##EQU11## Hence, by virtue of 7.6.31, the normal component of H_(e) satisfies the boundary condition: ##EQU12## at each point P of S_(h). The second term on the left hand side of equation 30 is the normal component of the intensity generated at P by the surface charge density σ. The symbol ρ denotes the distance of P from a point of S and the point Q whose charge m is located. As indicated in FIG. 3, the gradients of ρ⁻¹ are computed at point P. By virtue of equation 25, equation 30 transforms into the boundary equation ##EQU13## The integration of each term of equation 31 over the closed surface S_(h) yields: ##EQU14## where Ω)Q) is the solid angle of view of the closed surface S_(h) from the point Q where charge m is located. If point Q is outside of S₂, then

    Ω(Q)=0                                               (34)

Hence, by virtue of equations 32, 33 and 34, the integration of equation 61 over S_(h) yields: ##EQU15## which reflects the fact that the material of permeability μ immersed in the magnetic field generated by external sources is going to be polarized by the field, but it cannot acquire a non-zero magnetic charge.

In the limit μ=∞, S_(h) becomes an equipotential surface at a potential Φ_(h), whose value is determined by the solution of boundary equation 31. At each point P of S_(h), Φ_(h) is the sum of the potential generated by the charge distribution σ and by point charges m in a uniform medium of permeability μ₀. Thus, Φ_(h) must satisfy the equation: ##EQU16## where σ is given by the solution of equation 31. Since equation 35 is the direct consequence of equation 31, in the limit Φ→∞ the variables σ and Φ_(h) can be determined by the solution of the system of equations 35 and 36.

In the integral on the left hand side of equation 36, the distance ρ is zero for the element of charge σadS_(h) located at the point where the potential is computed. However, as long as σ is finite, the integral does not exhibit a singularity. Consider a circle on surface S_(h) of small radius and with center at P. For r→0, the potential due to the surface charge within the area πr² is ##EQU17##

A ferromagnetic material is characterized by a large value of its permeability. In the limit: ##EQU18## The normal component of H_(e) on the surface S_(h) may be written in the form: ##EQU19## where H_(eo) is the field intensity in the limit μ=∞ and factor G is a numerical factor that depends upon the geometry of S_(h). The G is a function of the position of the point P. By virtue of equations 29 and 30, the surface charge density σ(P) may be written in the form:

    σ(P)=σ.sub.∞ (P)+dσ                (40)

where σ.sub.∞ is the solution of equation 31 in the limit μ=∞. By virtue of equation 39,

    σ.sub.∞ =μ.sub.0 H.sub.e0                   (41)

Thus equation 40 yields: ##EQU20## By substituting the value of σ given by equation 40 in equation 31: ##EQU21## and by virtue of equation 42, function G satisfies the equation ##EQU22##

Once the value of dσ has been obtained by solving equation 43, the potential dμ generated inside surface S_(h) can be computed: ##EQU23## Thus, the magnetic induction B inside S_(h) is ##EQU24## i.e. in the limit of equation 38, the magnetic induction inside S_(h) is independent of μ and is determined only by the distribution of σ.sub.∞ and the geometry of S_(h).

In some particular case G is independent of the position of P, in which case dσ is proportional to σ.sub.∞, and the field generated by dσ, i.e. the external field in the absence of the medium of permeability μ.

As an example consider a cylinder of radius r₀ and permeability μ immersed in a uniform field of intensity H₀ perpendicular to the axis of the cylinder. Assume the polar coordinate system (r,Θ), where r is the distance from the axis of the cylinder and Θ is the angle between r and the direction of H₀. The radial component of the magnetic field is ##EQU25## and the surface charge density σ is ##EQU26## Thus in the limit (27)

    σ.sub.0 =2 μ.sub.0 H.sub.0 cosΘ             (49)

and ##EQU27## Thus the intensity δH of the field inside the ferromagnetic material is ##EQU28##

Numerical Solution

With the exception of some elementary geometries and distribution of magnetization like, for instance, a structure of concentric cylindrical or spherical layers of uniformly magnetized media and uniform materials, eqs. (12) and (13) cannot be solved in closed form, requiring numerical integration. This is accomplished by replacing in eqs. (12) and (13) the integrals with sums over small elements of surfaces of the ferromagnetic materials and the volume of the magnetized material. Thus, eqs. (12) and (13) transform to ##EQU29## where σ_(im) is the average value of the surface charge density in the element of surface δV_(n) · and (v·J)_(n) is the average value of the divergence of J in the element of volume δV_(n). The value ρ_(h),n is the distance between the center of an element of surface δS_(h) and the center of the element of volume δV_(m). The value ρ_(him) is the distance between the centers of elements of surface δS_(h) and δS_(im). The value Φ_(h) is the potential computed at the center of each element of surface δS_(h). Thus in the approximation of eqs. (39) and (40), the condition of constant potential is imposed only at a number of selected points equal to the number of surface elements. The potential is allowed to fluctuate between these points about the average values Φ_(h). The amplitude of the fluctuations decreases as the dimensions of the elements of the surface decrease.

As an example, apply eqs. (39) and (40) to the computation of the field in the two-dimensional problem of a strip of infinite magnetic permeability located in a uniform field as shown in FIG. 4, where the axis z coincides with the center of the strip. Assume that the uniform field is oriented in the positive direction of the axis y. If the potential is assumed to be zero on the plane y=0, the scalar potential of the uniform field is

    Φ=-H.sub.0 y,                                          (54)

where the positive constant H₀ is the intensity of the field. Because of symmetry, the potential of the strip must be equal to the value of the potential on the plane y=0, independent of the angle between the field and the plane of the strip. Thus in eq. (40)

    Φ.sub.h =0                                             (55)

The right hand side of eq. (40) corresponds to the potential at each point of the strip due to an external distribution of magnetization that generates the uniform field. Thus eq. (40) reduces to ##EQU30## where ρ is the distance of the m-th element of surface δS_(m) and a point P of the strip, and y is the ordinate of P.

The left hand side of eq. (43) can be readily integrated along the z coordinate. For a strip of infinite length, each element of surface of an infinitely long strip of infinitesimal width dζ generates a potential dΦ at a point P of the strip ##EQU31## where Φ is an arbitrary constant and r is the absolute value of the distance of P from the strip of width dζ:

    r=|ζ-τ|                         (58)

where ζ and τ are the distances of dζ and P from the center of the strip.

The numerical solution of eqs. (39) and (43) proceeds by dividing the width 2τ₀ of the strip in 2n equal intervals and by computing the left hand side of eq. (43) at the center of each interval. By virtue of eq. (28), if the number 2n of intervals is sufficiently large, one can neglect in each interval the contribution of the charges within the same interval.

Because of symmetry, the surface charge density satisfies the condition

    σ(-y)=-σ(y)                                    (59)

Thus eq. (39) is automatically satisfied and the values of σ(y) are the solutions of the system of n equations in the n variables σ_(m) ##EQU32## where coefficients a_(h),m are ##EQU33## for h≠m and ##EQU34## for h=m. In eqs. (47) σ_(m) is the average value of σ in the interval where the center has the coordinates ##EQU35## If α=π/2, i.e., if the external field is perpendicular to the strip, the solution of eq. (47) is

    σ.sub.m =0                                           (64)

for all values of m and no distortion of the field is generated by the strip. Thus the non zero value of σ_(m) is determined only by the field component parallel to the strip.

FIG. 5 shows the solution of the system of eqs. (47) for n=20. The plotting of the equipotential lines generated by the charge distribution of the strip is shown in FIG. 6. As expected, for Φ→0, the equipotential lines become circles that pass through the origin of the coordinates and with center located on the line ##EQU36## FIGS. 7 and 8 show the equipotential lines of the field around the strip in the two cases α=0 and α=π/4. In both cases the external equipotential lines Φ=0 intersect the strip at an angle π/2.

Once the field has been computed in the limit μ=∞, the field distortion generated by a small value of μ₀ /μ is obtained by the numerical solution of eq. (27). This is done by dividing S_(h) in a number n of small elements of surfaces δS_(m). Eq. (27) transforms to ##EQU37## where δσ so is the average value of δσ on the element of surface δS_(m), n_(k) is the unit vector perpendicular to the element of surface δS_(k), ∇_(k) is the gradient computed at a point infinitely close to the element of surface δS_(k) and inside S_(h), and ρ is the distance between the centers of δS_(k) and δS_(m). Thus eqs. (53) are the n equations in the n variables δσ_(m).

The system of eqs. (12) and (13) provides the exact solution of the field generated by an arbitrary distribution of remanences in a transparent medium (χ_(m) =0) limited by a number of surfaces of infinite magnetic permeability materials and arbitrary geometries.

In a structure of media of uniform values of χ_(m) and μ, the solution of eqs. (23) and (24) is proportional to χ_(m) and the solution of eq. (32) is proportional to μ₀ /μ. Thus the scalar potential at each point P of the magnetic structure is ##EQU38## where Φ₀ is the potential in the ideal case χ_(m) =0 and μ₀ =0, and ψ₁, ψ₂ are functions of position which are determined by Φ₀, independent of χ_(m) and μ₀ /μ. Usually, the rare earth magnetic materials exhibit values of the order of 10⁻² and the linear range of the characteristic values of μ₀ of the order of 10⁻³ or smaller.

Thus, outside of the ferromagnetic components of the structure one can expect the demagnetization characteristic to be the dominant factor in the field perturbation.

An example of the numerical solution is the field computation in the two-dimensional problems of a high permeability material whose cross section is the equilateral hexadecagon shown in FIG. 9 with sides tangent to an ellipse with 2:1 ratio between axes. The external uniform field of intensity H₀ is oriented at an angle π/4 with respect to the axis of the ellipse. The equipotential surface Φ=0 of the external field is assumed to contain the axes of the polyhedron.

The field corresponding to a finite (μ₀ /μ=0.5) magnetic permeability, computed according to equation (45), is plotted in FIG. 10.

An example of multiplicity of high permeability components is the two-dimensional structure shown in FIG. 11. The two lined rectangular areas represent the magnetic material uniformly magnetized in the direction of the y axis. The heavy lines represent the cross-sections of four components of zero thickness and infinite permeability.

The field configuration derived from the numerical solution of equation (31) is shown in FIG. 12. In this figure the equipotential lines are plotted in the first quadrant of the structure of FIG. 11. The numerical solution is shown for y₁ =2y₀ =x₀. The x axis is a Φ=0 equipotential line within the region of the magnetized material that intersects the x axis at a point X that becomes a saddle point of the equipotential lines. The numerical values of the potentials are Φ₁ =-Φ₂ =-0.248, Φ₃ =-Φ₄ =0.277.

FIG. 13 illustrates the field configuration in the case of separation of the inclined sides. As can be seen, the surfaces acquire a potential different from the configuration shown in the previous example.

If S₃ and S₄ are assumed to be connected to each other at infinity, FIG. 11 may be considered as the ideal schematization of a yoked magnet. In this case both Φ₃ and Φ₄ are zero. FIG. 14 shows the equipotential lines of the field computed within the structure and FIG. 15 shows the field outside. Point Y on the y axis is a saddle point of the field configuration. The field in the region between surfaces S₁ and S₂ has approximately the same magnitude as the field within the magnetized material. This is the result of enclosing the magnetized material within the yoke formed by the surfaces S₃ and S₄.

FIGS. 16, 17 and 18 are self explanatory flow diagrams illustrating an example of the invention. As noted, FIG. 17 constitutes a continuation of FIG. 16, and FIG. 18 constitutes a continuation of FIG. 17.

While the invention has been disclosed and described with reference to a single embodiment, it will be apparent that variations and modification may be made therein, and it is therefore intended in the following claims to cover each such variation and modification as falls within the true spirit and scope of the invention. 

What is claimed is:
 1. A method for constructing a permanent magnetic structure with linear magnetic characteristics, comprising specifying dimensional parameters of a permanent magnetic structure having interfaces between magnetized regions, predetermined remanence and susceptibility characteristics, determining the surface charges σ at each interface of the magnetized regions, dividing the surface of the structure into a plurality of predetermined surface regions with each of said regions having a defined point, determining the distribution of said surface charges on all of the interfaces, computing the surface charges σ, then computing the field everywhere using the calculated surface charges, then repeating said steps of specifying dimensional parameters, determining surfaces charges, dividing, and determined the distribution of said surface charges until said computed field is a determined value, and then fabricating a permanent magnetic structure in accordance with the last specified dimensional parameters.
 2. A method for constructing a permanent magnetic structure comprised of components of both magnetic and ferromagnetic materials, with linear magnetic characteristics, comprising specifying dimensional parameters of a permanent magnetic structure having interfaces between magnetized regions, assuming infinite permeability of the ferromagnetic components, determining surface charges at each said interface, formulating a set of linear equations of said structure in terms of the scaler potential, determining charge elements of said structure from said charge equations, determining the field of said structure from said elements, then repeating said steps of specifying dimensional parameters, determining surface charges, formulating a set of linear equations, determining charge elements, and determining the field until the determined field is a desired value, and the fabricating said permanent magnetic structure in accordance with the last specified dimensional parameters.
 3. The method of claim 2 wherein said step of determining the field of said structure comprises directly determining the expansion of the magnetostatic potential.
 4. A method for constructing a permanent magnetic structure comprised of components of both magnetic and ferromagnetic materials, with linear magnetic characteristics, comprising specifying dimensional parameters of a permanent magnetic structure having interfaces between magnetized regions, assuming finite permeability of the ferromagnetic components, determining the surface charges at each said interface, dividing the interfaces into a plurality of surface regions, formulating a set of linear equations expressing the surface charge elements of said regions in terms of scaler potentials, determining charge elements of said structure from said equations, determining the field of said structure from said charge elements, then repeating said steps of specifying dimensional parameters, determining the surface charges, dividing the interfaces, formulating a set of linear equations, determining charge elements and determining the field of the structure until said determined field has a desired value, and then fabricating a permanent magnetic structure in accordance with the last specified dimensional parameters.
 5. The method of claim 4 wherein said step of determining the field of said structure comprises directly determining the expansion of the magnetostatic potential.
 6. A method for constructing a permanent magnetic structure comprised of components of both magnetic and ferromagnetic materials, with linear magnetic characteristics, comprising specifying dimensional parameters of a permanent magnetic structure having interfaces between magnetized regions, assuming finite permeability of the ferromagnetic components, determining the surface charges at each said interface, dividing the interfaces into a plurality of surface regions, formulating a set of linear equations expressing surface charges of said structure in terms of the vector field intensities, determining unknown charge elements of said structure from said equations, determining the field of said structure from said charge elements, then repeating said steps of specifying dimensional parameters, determining the surface charges, dividing the interfaces, formulating a set of linear equations, determining unknown charge elements, and determining the field, until a predetermined field is determined, and then fabricating a permanent magnetic structure in accordance with the last specified dimensional parameters.
 7. The method of claim 6 wherein said step of determining the field of said structure comprises directly determining the expansion of the magnetostatic potential. 